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Abstract 

(N 

K*" I We present analytical and numerical results for static, spherically symmetric solu- 

' tions of the Einstein- Yang-Mills-Higgs equations corresponding to magnetic mono- 

. poles and non-abelian magnetically charged black holes. In the limit of infinite Higgs 

mass we give an existence proof for these solutions. The stability of the abelian ex- 
■ tremal Reissner-Nordstr0m black holes is reanalyzed. 
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In the present paper we continue our investigation of magnetically charged nonabelian black 
holes, resp. gravitating magnetic monopoles, started in Q, referred to as I. The latter are 
I solutions of an Einstein- Yang-Mills-Higgs (EYMH) theory with the Higgs field in the adjoint 

representation. In the following we restrict ourselves to an SU(2) gauge theory as in I. 

Gravitating magnetic monopoles were first considered in Q and later studied numerically 
together with black holes in some detail in |]|, ^, Q|. The numerical studies revealed a number of 
interesting phenomena for both of the cases when the strength of the gravitational force is either 
very weak or becomes comparable to that of the YM interaction. The relevant dimensionless 
parameter characterizing the relative strength of the gravitational interaction is the mass ratio 
a = Mw/fl'-/Vfpianck where Mpianck = ^/VG is the the Planck mass and Mw is the mass of 
the YM field (due to the Higgs effect). It has been found that regular monopoles exist only for 
a < Omax of order one. The precise value of Omax depends on the second dimensionless mass ratio 
(3 = M}i/M-\!\i where Mh is the mass of the Higgs field. In addition to the gravitating monopole 
solution tending to its flat limit for a — > a discrete family of radially excited monopole solutions 
was found for arbitrarily small a (and a < \/3/2) Q having no counterpart in flat space. 

Another new feature without analogue in flat space is the existence of nonabelian magneti- 
cally charged black holes. They were found to exist in a certain bounded domain in the (a, r^) 
plane where denotes the radius of the black hole horizon. The precise form of this domain 
again depends on (3. 

^Permanent address: Research Institute for Particle and Nuclear Physics H-1525 Budapest 114, P.O.Box 49 



1 



Whereas the emphasis in I was mainly on numerical results, the present work is devoted 
more to to analytical considerations and a qualitative understanding of the numerical results. 

It would be clearly desirable to give analytical proofs for the existence of the monopole 
solutions found by numerical integration as it was achieved |5|, for the discrete family of 
solutions of the EYM theory found by Bartnik and McKinnon [0] and the related black holes 
[^. One of the main results of this paper is an existence proof of magnetically charged non 
abelian black holes, resp. magnetic monopoles in the limiting case /3 — > oo when the Higgs 
field is frozen to its vacuum expectation value. In this case the analysis of the problem is 
considerably simplified while many features of the full problem (i.e., the one with a dynamical 
Higgs field) do carry over. Although the corresponding spherically symmetric ansatz for the 
Higgs field is singular at the origin, there exist still 'quasi regular' solutions with finite mass. 
The corresponding gravitational field develops a 'conical singularity' at the origin, whereas the 
YM field is still differentiable. The angular deficit leading to this conical singularity is the same 
as the one found for 'global monopoles' 1^]. The only difference is that in the latter case it has 
been transferred to large values of r using I/Mh as length scale instead of 1/Mw used here. 

The maximal value of the mass a. (measured in units of the Planck mass) for which the quasi 
regular solution exists was found to be l/\/2 This is exactly the value where the angular 
deficit becomes 47r, i.e., the conical singularity changes its nature. The behaviour of the solution 
at the origin is then rather different from the one for smaller values of a. A similar problem 
arises, when one considers black holes with a degenerate horizon. In fact, the quasi regular 
solution emerges in the limit r/^ ^ of such extremal black holes. Even for local solutions the 
proof of existence for extremal black holes (and their numerical construction) requires certain 
modifications of the methods used in I. These will be described in Sect. |^, after discussing the 
field equations in Sect. ^ and the necessary boundary conditions in Sect. ^ 

Sect. ^ is devoted to an analysis of the global behaviour of all solutions with the boundary 
conditions at r = 0, resp. discussed in Sect. ^ Methods and results of this analysis, leading 
to a complete classification, follow closely those of @. As a consequence of this classification the 
existence of global black hole solutions for the massive EYM theory (/3 = oo) is demonstrated 
for a certain domain of the (a, r/j) plane. Analyzing the limiting behaviour of the solutions when 
the horizon degenerates, resp. tends to zero the existence proof can be extended to extremal 
black holes, resp. quasi regular solutions. 

In Sect. |6|we discuss the problems involved in the numerical solution of the relevant boundary 
value problem with particular emphasis on the difficulties arising from the singular nature of 
the boundary points. Furthermore we try to give a qualitative understanding of the various 
mechanisms delimiting the existence region of global (regular and black hole) solutions as a 
function of the parameters a and /?. 

Some numerical solutions are also presented to illustrate the usefulness of the techniques 
described in Sect. ^ in a case, which could not be dealt with previously (/3 = oo, a = l/\/2). 

In the last section we address the stability of the extremal Reissner-Nordstr0m solution. We 
demonstrate that it becomes unstable (independently of /?) for precisely a = \/3/2 and not for 
a = 1 as it has been claimed previously |Q, |T^. This result complies well with the numerically 
observed bifurcation of the nonabelian black holes of radius = a with the extremal RN 
solution for all values of /?. 



2 Ansatz and Field Equations 

We are interested in static, spherically symmetric solutions of the EYMH equations. In this case 
the metric tensor of the space-time can be parametrized as [12| 



(1) 
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where dO,'^ = d0^+sin^ Odip"^. Whereas the functions e'^, resp. 47rr^ have a geometrical significance 
as the length of the time-translation Killing vector, resp. the surface area of the invariant 2- 
spheres the function is gauge dependent, i.e., depends on the choice of the radial coordinate 
R. A very convenient choice is obtained putting = r/R corresponding to isotropic coordinates 
for the 3-spaces t =const. 

For the SU(2) Yang-Mills field we use the standard minimal spherically symmetric 
(purely 'magnetic') ansatz 

W^Tadx^" = W{R) {Tide + T2 sin Odip) + cos ed^p , (2) 

where Ta denote the generators of SU{2) and for the Higgs field we assume the form 

$» = H{R)n'' , (3) 

where is the unit vector in the radial direction. The reduced EYMH action can be expressed 

as 

S = - j di?e(^+^) (1 + e-''\{r'f + v\r'')') - Wi - V2] , (4) 



with 



and 



V.J^^liHr. (5) 



r 



V2 = + ^{H'- a'f + W'H' . (6) 

Through a suitable rescaling wc have achieved that the action depends only on the dimcnsionless 
parameters a and (3 representing the mass ratios a = M\\/\fGlg = M^y/^'Mpi and /3 = Mh/Mw 
(g denoting the gauge coupling). 

In order to derive the gravitational field equations we have to employ independent variations 
of the functions v, A and r, inserting the gauge choice = r/R afterwards. 

The resulting equations are 



i (1 - {—f - 2R'u''-) + R'V, - = , (7a) 
1 (1 _ (:?^)2 _ 2-{Rt')') - R'Vi -V2 = 0, (7b) 



R{ ) -h-Re {Rev)^Rr— hr— — = 0, (7c) 

r or or 

T> \Al^ 1 

Ri^e^'W')' -e'rW^ +H^)=0, (7d) 

RiRre^H')' - e''rH{2W + ^{H^ - a^)) = . (7e) 

In order to obtain an autonomous first order system we introduce the coordinate r = In i? (with 
'= d/dr) and new variables for the first derivatives through 

N=- , K = iy + N , U= — , and V = H . (8) 
r r 

Thus we obtain the system of equations 

r = rN , (9a) 

N={k- N)N -2U^ -V^ , (9b) 
g2 2 

k=l- + 2U^ - ^-^{H^ - a^f - 2H^W^ , (9c) 
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W = rU, (9d) 

jj^WiW^^2}l+rH^W-{K-N)U , (9e) 
r 

H = V , (9f ) 

a2 2 

V = ^-^{H'^-a^)H + 2W'^H-KV , (9g) 
together with the constraint 

2kN = 1 + iV^ + 2C/^ + y2 - 2V2 , (10) 

derived from Eq. (0a). 

A certain simplification occurs taking the formal limit /? — > 00 of infinite Higgs mass. In this 
limit the Higgs field is frozen at its 'vacuum' value H = a and at the same time the last two 
of the field equations are discarded. As seen from the ansatz (^) the Higgs field is singular at 
r = 0, but this singularity is not directly visible in the radial field equations 

r = rN , (11a) 

N={k- N)N -2U'^ , (lib) 

k = l- + 2U'^ -2a'^W'^ , (11c) 

W = rU, (lid) 

U=^^^^^^^—^ + ra^W-{K-N)U, (lie) 



and the constraint 



(I — w2\2 

2kN = 1 + N'^ + 2U'^ - ^ - 2a'^W'^ . (12) 



These equations may be interpreted as describing a massive YM field (of mass a). For a = 
the equations coincide with the Eqs. (49,50) of Q. 

Special solutions of Eqs. @, resp. (|Tl| ) are the abelian Reissner-Nordstr0m (RN) black holes 
with (r/i denoting the radius of the horizon) 



W = , H = a , = 1- — ^ + — . (13) 



rrfi r 

Expressing r as a function of r yields 



,9''" 1 • , 9 ,r (r? — l)sinhr , 

r = rh cosh^ smh^ - , N = ^ , k = coth r , (14) 

2 rh 2 ' 2rrh ^ ' 

for the regular RN solution with > 1 and 

r = l + e^ N = -^-^, K = l, (15) 
1 + e ^ 

for the extremal RN solution with r/^ = 1. 



3 Boundary conditions 

The field Eqs. (^ have to be complemented with suitable boundary conditions. For physical 
reasons we are only interested in solutions with a regular origin r = or black holes. In both 
cases the corresponding boundary points are singular points of the field equations. Whereas for 
r = this is obvious from Eqs. (^ it is less obvious in the case of black holes. Using Eqs. (^) 
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and = r/R one obtains the relation N = e~^dr / dR. Since has to vanish on the horizon 
of the black hole also N vanishes there. This, however, leads to a divergence of k in view of the 
constraint Eq. ([lO|). 

For asymptotically flat solutions we also have to impose boundary conditions at r = oo, 
another singular point of the equations. 

Usually the regular solutions at some singular point are determined by a number of free pa- 
rameters strictly smaller than that for a non-singular point. A standard procedure to determine 
these free parameters is to expand the solution in a formal Taylor series in the independent vari- 
able and read off the undetermined coefficients. If this Taylor series converges it defines a family 
of local solutions parametrized directly at the singular point. There are, however, cases where 
this method fails, e.g., for black holes with a degenerate horizon. Yet, there is a more general 
method to obtain local solutions at singular (fixed) points of 'dynamical systems' (i.e., systems 
of first order ODE's). One tries to find a linearization of the differential operator which allows 
to characterize the behaviour of solutions near the singular point. In the case of a hyperbolic 
fixed point regular solutions of the linearized system lie on a linear subspace ('stable manifold'). 
Regular solutions of the full non-linear system lie on a corresponding curved stable manifold of 
the fixed point, whose existence is guaranteed by standard textbook theorems (e.g., [|l3|, Q). 
The method to determine the stable manifold through a system of integral equations solvable 
by iteration yields a viable numerical method to obtain local solutions even in cases, where the 
Taylor series approach fails. 

After these general remarks we turn to a detailed discussion of the various boundary condi- 
tions. The behaviour of regular solutions at r = was already treated in I, although with the 
coordinate choice r = R (Schwarzschild coordinates). One finds the asymptotic behaviour 

W{r) = 1 - + O(r^) , (16a) 
H{r) = ar + 0{r^) , (16b) 

N{r) = 1 - l(a2 + 462 + + 0{r^) , (16c) 

where a and h are free parameters. In view of the preceding discussion we interprete the ex- 
pansion of W etc. in terms of r as an approximation to the 'stable manifold', i.e., as a relation 
between dependent variables. In this case the parameters a and b are defined as the limits of 
H/r and (1 — W'^)/r'^ for r — > — oo, i.e., at the singular point. 

Next we turn to the case (3 = oo with the field Eqs. (P). Also this case has already been 
discussed in I, where we found 

2 

W{r) = l-^— + ... + br^ + o(rT) , (17a) 

1 - 3 



N{r) = Vl - 2a2 1 + — ^ _ + . . . + br^ + o(r^) , (17b) 



the dots representing even polynomials in r of degree smaller than 7. As long as < 1/2 
the stable manifold can still be parametrized by b, although to define it through the limit of 
(1 — W'^)/r'^ one has first to substract a polynomial of degree smaller than 7. We remark that 



whenever 7 is an even integer larger than 2, i.e., ca = \J \ ~ l/2m(2m — 1) with m = 2,3, .. . 
logarithms appear in Eqs. (p!7|). 

The most important new feature is the behaviour of at r = 0. Whereas A^(0) = 1 for all 
finite /?, now A(0) = Vl - 2q^ implying an angular deficit and a curvature singularity at the 
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Figure 1: The curve iVi = for /? = oo and = 0, 0.3, 0.5, 0.6, 0.75, 1, and 1.5 
Figure 2: Some regions of the (r, W) plane for a = 0.75 and (3 = co 



origin. It is induced by the singularity of the Higgs field at r = 0. The asymptotic expression for 
N{r) shows that a has to be restricted to < a < l/\/2- In the limit a l/\/2 the expression 
for 7 and the conditions (17) become meaningless. This is a case where the 'stable manifold' 
itself is well-defined, however, its parametrization at the singular point breaks down. A precise 
formulation of this case will be given in the next section. 

The boundary conditions for regular black holes at the horizon were already given in I in 
Schwarzschild coordinates. For later convenience we repeat them here using the coordinate r. 
As already mentioned the function k{t) is singular at the horizon (which we choose to put at 
r = 0), but k{t) — 1/t turns out to be regular. One obtains 



r{T)=rh \l + Ni 



where 



iV-i 



N{t) 
W{t) 
H{t) 

_ 1 
~ 2 ' 



_2 



Wh + rhWi— + 0{T' 



2\h 



W^ 



Th dV2 

4 dW 



1 m 

2 dH 



(18a) 
(18b) 
(18c) 
(18d) 

(19) 



and Vh, Wh, and Hh are undetermined parameters. For (3 = oo one has to substitute H = a, 
Hi = and omit Eq. (^d). In order to have a regular horizon we have to require A'^i to be 
positive. For given r^, a, (3 the inequality A'^i > determines some bounded (allowed) region in 
the iWhiHfi) plane. In the limiting case (3 = oo the expansions simplify in an obvious way and 
the allowed region of the (r/j, Wh) plane is bounded by the quartic curve A^i = 0. In Fig. |l| these 
quartic curves are shown for various values of a, the allowed region is shown as C = C+ U C_ 
for a = 0.75 in Fig. || 

If /? > 2 the regions in the {Wh, Hh) plane allowed by the condition iVi > shrink to points 
for certain values of a and rh- In this case Wh and Hh are fixed by the conditions 



{Wl - If 



+ 2WlHl + '^rl{Hl 



4 



a 



2^2 . 



^ + rlHl. 



^1 , 
0, 



(20a) 
(20b) 
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Figure 3: The curves a^rh vs. for /? = 3, 4, 6, 10, 20, and oo 



trl{Hl-a') + Wl = 0. (20c) 

Nevertheless there still exists a two parameter family of local solutions, but members of this 
family cannot be characterized by their behaviour at the horizon. The horizon is degenerate 
since Eqs. (|o|) imply Ni = Wi = Hi = 0. 

It is convenient to give the solution of Eqs. (20) in parametric form: 

rlHl = l-Wi, (21a) 



a^rl = l-il-^)Wl, (21b) 
4 



d = 'i--a-^)W^. (21c) 

When Wh increases from to 1 for some fixed (3 > 2, then Hfi decreases from 1 to 0, 
decreases from 1 to 2//3, and a first decreases from 1 to a minimal value and then increases to 1; 
the minimal value Omin = 

V(/3 + 2)/2/? is reached for Wh = ^/(3/{(3 + 2) where rh = 2/V^T2 
and Hh = 1/V2. Solving Eqs. (|2l|) for a in terms of rh one obtains a function ae{rh) given by 



1 



^(l_4//?2)(l_,2) 



al = L_ __ 1_ . (22) 

This function is shown in Fig. ^ for various values of /?. 

The conditions A^i = Wi = Hi = for a degenerate horizon have some more solutions 
besides Eqs. (|^), discussed in more detail in Sect. ^. 

Finally let us consider the b.c. when r oo. Imposing 

iV(r) = l-^ + 0(l) , (23) 
ensures asymptotic flatness. Then from the field eqs. one easily finds: 

Wir) = Ce""'^ (^1 + O (^i^ ^ , (24a) 

H{r) = a- B^-^ (^1 + (^^^^ , fov P <2 , (24b) 

^M = «-^;32 r^^^i^ + Oi-)) , fov I3>2 , (24c) 



a(/3^ — 4) V V?' 
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where a = r + M Inr and M, B, and C are free parameters, i.e., there exists a three-parameter 
family of asymptotically flat solutions. The asymptotic behaviour of the Higgs field H{r) depends 
crucially on /3. When /? > 2 its asymptotic decay is dominated by the non-linear terms in 
Eqs. (^) and it is not possible to fully parametrize the stable manifold at the singular point 
r = oo. 



4 Local Existence 

In the preceding section we have analyzed the behaviour of solutions of Eqs. (p|,p!o[) or (11,12) 
that remain regular at one of the singular points (r = 0, r = oo, or a horizon). We have 
found that there are families of such solutions depending on one or several free parameters. The 
arguments were, however, based on a formal power series expansion. In order to show that such 
solutions exist as real analytic functions we apply standard techniques to linearize non-linear 
systems at a singular point. 

For regular black holes k has a pole at the horizon for finite r and in Eqs. (^) we have chosen 
T = at the horizon. In all other cases discussed in Sect. ^ the dependent variables have a limit 
as r — > iboo such that all the right hand sides of Eqs. (P), resp. (11) vanish. The behaviour of 



solutions near such a singular or fixed point can be characterized by a 'stable manifold' (see, e.g., 
|l3|| p.330ff). Consider a system of autonomous first order differential eqs. for m + n functions 

y = (Ui, . . . ,Um,Vl, ... ,Vn) 

^=Au + fiy), ^ = -Bv + g{y), (25) 

where all eigenvalues //, resp. A of the constant matrices A, resp. B have a positive real part, and 
f{y)j 9{y) cire 0(||y|p) and analytic for ||y|| sufficiently small. There exists an n-dimensional 
stable manifold (for t — > +oo) of initial data y(0) = {ip{v),v) such that y ^ as t — > -|-oo. 
Identifying solutions with the same orbit (only differing by a shift in t) there exists an (n — 1)- 
parameter family of such equivalence classes of solutions on the stable manifold. The functions 
^p^{v) are 0(||'y|P) and analytic for sufficiently small. For such v and t > the solutions are 
analytic in v and t. Due to the presence of non-linear terms this result does not, however, imply 
that e^^v{t) has a limit for t oo. For solutions off the stable manifold the distance ||n — ■i/;(f)|| 
increases exponentially as long as ||y|| is small. 

Similarly there exists an m-dimensional stable manifold (for t — > — oo) of initial data y(0) = 
{u, ip{u)) such that y — > as t — > — oo. 

It is sometimes possible to characterize solutions on the stable manifold by their behaviour 
at the singular point. In we have derived the following result for these cases (see also 
|14| p.304ff): Consider a system of first order differential eqs. for m + n functions y = {u,v) 



s^ = sfi{s,y), i = l,...,m, (26a) 
as 

dv ' 

s—^ = -XiVi + .sgi{s,y), i = l,...,n, (26b) 
as 

with constants Aj > and let C be an open subset of such that /j, gi are analytic in a 
neighbourhood of s = 0, y = (c, 0) for all c £ C. There exists an m-parameter family of local 
solutions yds) analytic in c and s for c G C, \s\ < so{c) such that yc(0) = (c, 0). 

We have used this result in [|l| to show that there exists a two-parameter family of solutions 
of Eqs. (PJlOj) with regular origin where W, H, and N are analytic in r, a, and b and satisfy 
the boundary conditions (|l^). Finally r is an analytic function of r due to Eq. (|9|a) and k as 
computed from the constraint (|To| ) is analytic as well. 
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Let us next consider the 'quasi regular' solutions of Eqs. (11,12) with a conical singularity 
at the origin. Assuming a < 1/ \/2 we define v = Vl — lo?-, introduce 

{W -\) r 
r 4 



W 



V 1 



N = N 



and use the constraint ( |l^ ) to express k in terms of r, W, U, and A^. Taking the linear 
combinations X+ = {'y — 1)W + U and X_ = 'yW — U with 7 from Eq. (^c) we obtain a 
system of equations in the form (^) for t = t, u = (r, v = {N,X^ ) with eigenvalues 

/u = (z^, (7 — 1)1^), A = [v^^v). Hence there exists a two-dimensional stable manifold for r — > 
— 00, i.e., a one-parameter family of quasi regular solutions that can be characterized, e.g., by 
X+(r) = —b{r)A'^~^\ As long as 7 < 4 the system can be written in the form (|2^), and hence 
b{r) has a limit for r ^ 0; for larger 7 higher powers of r must be subtracted from X± in order 
to parametrize the solutions at the singular point r = (compare Eq. (p!7|a)). 

To show the existence of solutions satisfying the black hole boundary conditions ( [l8|) we 
observe that Eqs. (^) can be written in the form ( p^ ) with s = t, u = {r, W, H, 2kN — 2C/^ — V"^), 
V = {U, V,K — and A = (1, 1, 2). One free parameter is fixed by the constraint (p^), hence 
there exists a three-parameter family of local solutions analytic in rh, Wh, H^, and r. These 
solutions describe a regular horizon as long as {rh, Wh, Hh) are chosen such that A"i > 0. Since 
W , H, and A^^ depend only on and A^i 7^ they may be expressed as analytic funtions of 
r. Similarly there exists a two-parameter family of local solutions of Eqs. (|ll,12) analytic in r^, 
Wh, and r (or r). 

A degenerate horizon occurs if (r, A^, k, W, U, H, V) tend to a fixed point of Eqs. (^) with 
N = U = V = OasT^ —00. For such a fixed point A^i, Wi, and Hi in Eq. (|l9| ) have to vanish, 
furthermore = —rdV2/dr requiring rdV2/dr < 0. (For rdV2/dr > there are infinitely many 
zeros and poles of k instead of a fixed point.) There are four possible types of fixed points: 



1. rh = l,Wh = 0, and Hh 

2. solutions of Eqs 

3. rh = Wh = 

4. l-rl + a^p^rfl/4 



a. 



1, and Hh = 0, 
= 0,Wh = 0, and Hh 



0. 



In order to analyze the behaviour near the fixed points we introduce r = r — rh, k = k — Kh, 
W = W — Wh, and H = H — Hh- Using the constraint (^) to express f in terms of the other 



quantities we obtain r = Nrh/Kh up to non-linear terms, 
have the form y = My with y^ = {N, W, U, H, V, k) and 



The resulting linearized equations 



/ 



M 







rh 










^WhHh 



\ 





rh 

-Kh 











2rhWhHh 
'-^ + 2W, 



h 








1 

-Kh 





\ 







-2Kh ) 



(27) 



where • denotes some non-vanishing matrix elements whose precise form is not required. The 
matrix M has a positive eigenvalue Kh and a negative one —2Kh. In order to determine the 
remaining four eigenvalues 7^ let us study each of the possible fixed points in some detail. 
First there is, for all a and [3, the fixed point r = K = l, W = 0, H = a with eigenvalues 



71,2 



13,4 ■■ 



lb iuj with 



(28a) 
(28b) 
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The two eigenvalues 73^4 due to the Higgs field are real and have opposite sign for a/3 7^ 0. The 
two eigenvalues 71^2 due to the Yang-Mills field are complex conjugate and have a negative real 
part for a < \/3/2, are both negative for \/3/2 < a < 1, and have opposite sign for a > 1. 

For a > 1 there exists a two-parameter family of local solutions with a degenerate horizon 
approaching the fixed point as r —00. For a < 1 there is a one-parameter family with 
W = 0. The extremal RN solution with W = and H = a is a member of these families. In 
addition for a > 1, resp. a < 1 there exists a two-, resp. three-parameter family of local solutions 
approaching the fixed point as r — > +00. For a < y/3/2 the asymptotic behaviour of W and U 
is dominated by the linearized equations, i.e., 

W ^ Ce-^/^ sm{ujT + 6) , (29) 

whereas W monotonously decreases for a > V3/2 and r sufficiently large. The two-parameter 
family of solutions for q > 1 is characterized by e^^'^W ^ as r ^ +00, whereas e^/'^W is in 
general unbounded for the three-parameter family of solutions for \/3/2 < a < 1. 

Next there are the solutions of Eqs (^), which turn out to be of interest only for j3 > 2. As 
long as rdV2/dr < there are two positive and two negative eigenvalues 



7i = ± 




+ K + '^rlHl ± W - ^rlHl + mlrlHl , (30) 



2 ^ 4 2 

i.e., the stable manifolds for r itcxD are both three dimensional. For given a and (5 and each 
set of values {rh,Wh, Hh) such that Eqs ( pO| ) are satisfied there exists a two-parameter family 
of local solutions with a degenerate horizon approaching the fixed point as r — > —00, and there 
exists a two parameter family of local solutions approaching the fixed point as r ^ +00. 

For W = 1 and H = there is no fixed point since rdV2/dr > 0. Finally, for W = H = 
and as long as rdV2/dr < all four eigenvalues 7^ have a negative real part. There exists only 
one degenerate black hole solution, the extremal RN solution in a de Sitter background with 
W = H = 0, and k — > —00, — > for some finite r. 

The degenerate horizons for the theory with (3 = 00 are similar to those for finite /?. Dis- 
carding the two rows and columns corresponding to H and V and substituting H = a in the 
matrix M (Eq. (^)) we obtain four eigenvalues Kh, —2Kh, and 71^2- There are, however, just 
two possible fixed points. 

For the fixed point r = k = 1, 1^ = the eigenvalues 71^2 are given by Eq. (p8|a) above. For 
a > 1 there exists a one-parameter family of local solutions with a degenerate horizon, for a < 1 
there exists only the extremal RN solution with W = 0. In addition for a > 1, resp. a < 1 
there exists a one-, resp. two-parameter family of local solutions approaching the fixed point as 
r — > +00 and again W oscillates with an exponentially decreasing amplitude for a < \/3/2. 

For the fixed point = 1 — W^, oP'r^ = (1 — M^^), k = cP'r there is one positive and one 
negative eigenvalue 

li,2 = -^±^^ + 2Wi. (31) 

Therefore there exists a one-parameter family of local solutions with a degenerate horizon and 
a one-parameter family of local solutions approaching the fixed point as r ^ +00. 

The general results about the existence of a stable manifold require that none of the eigen- 
values have a vanishing real part. This is no longer the case for a degenerate horizon with 
Kh = 0, e.g., at the point on the curve (^T|) where a has its minimal value. We are interested in 
particular in the case (3 = 00 where this point r/j = 0, Wh = 1, a = \/^/2 can also be seen as 
limiting case of the quasi regular solutions. Guided by a Taylor series expansion we introduce 

W-l U 1 r - y/VlN 1 _ \/3k 1 

^± = ~ ^ + 7 ^ -fKF ' ^ = 2 ' = — ' 

V2r 4 v96 j,^ j, ^ 
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and use the constraint (^) to express k in terms of X±, N and r. The resulting differential eqs. 
are (with ' = d/dr) 



X4= (--±^W± + /±, N' = --N + fo. (32) 

The 'non-linear terms' f±fi as well as R/r are bounded as long as r is small and X±/r, 
N /r are bounded. In order to show that there exists a one-parameter family of solutions of the 
non-linear system (|3^ ) we closely follow the treatment for the standard case (|25| ) in p.330ff. 
The linear system obtained by substituting f±fl = in Eqs. (]3^) has the solutions 



^+(0)(0 = (7)^"^^'"^^^+(^o) , X_(o)(r) = iV(o)(r) = , (33) 

with X-|_(ro) as a free parameter. A solution of the non-linear system with the same boundary 
conditions X+(ro) and ^-(0) = -/V(O) = satisfies the integral equations 

X^{r)=X^^o){r)-f (^)%-^(^-T)/+(t)(it , (34a) 



^-^"^^l ''^'^'''^f-(*)dt, (34b) 

N{r) = f^(^-^^ mdt- (34c) 

Iterating these integral eqs. with the solution (33) of the linearized system as zeroth approx- 
imation one can easily estimate that the successive approximation converge to a solution of 
Eqs. (^) provided ro and X+(ro) are sufficiently small. Apart from showing the existence 
of a one-parameter family of local solutions this method also allows to compute the solutions 
numerically by iterating a suitably discretized version of the integral Eqs. 



Let us finally consider the limit r — > oo with the two-parameter family (^,24 a) of asymp- 
totically flat solutions of Eq. (pT|). Here we introduce 



m(r) = -(l-iV^) , s(r) = r + m(r)lnr , Y±{r) = e'^'^^'^W , 
2 \aj 

where m is monotonously increasing and M = m(oo) is the total mass. Expressing k in terms 
of m, Y±, and r we obtain the differential eqs. (again with ' = d/dr) 



m 



' = fo, Y; = U, r = 2a.s'y_ + /_ , (35) 



with 'non-linear terms' /o, resp. f± decreasing at least as r~^, resp. r^^lnr as long as m and 
Y± are bounded. Together with the boundary conditions m(oo) = M, y+(oo) = C, y_(oo) = 
we obtain the integral equations 



m{r) = M-J fo{t)dt , (36a) 

U{t)dt , (36b) 

POO 

Y_{r) = - e-2"(^W-^(''))/-(t)'^i • (36c) 



Iteration of these integral eqs. with m(o) = M , l+(o) = and ^-(o) = as zeroth approximation 
yields a convergent series for any given M and C and sufficiently large r. Note that in Eqs. (|35| , |36|) 
we have used s(r) = r + m(r) In r and not cj(r) = r + Mlnr as in Eqs. (p^). This allows to 
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rewrite the integral eqs. in a form more suitable for our numerical procedure: Given m(ro) and 



Y'+(ro) for some sufficiently large value rg, a slightly modified version of Eqs. (36) allows to 
compute y_(ro) as well as the asymptotic parameters M and C. 

The method described above can be easily adapted to the three-parameter family of asymp- 
totically fiat solutions with finite (3 as long as /3 < 2. The situation is, however, slightly different 
for /? > 2, where the asymptotic behaviour of the Higgs field is dominated by non-linear terms. 
There still exists a three-parameter family of solutions but at least one of the parameters must 
be determined away from the singular point r = oo by an integral eq. of the type (| 



5 Global Existence for /3 = oo 

In this section we shall discuss the global behaviour of solutions of the theory with (3 = oo 



(frozen Higgs field) whose field equations (11) were derived in Sect. We restrict our discussion 
to solutions with either 'quasi regular' boundary conditions at r = (requiring a < or 
those for black holes. These solutions depend on one, resp. two parameters, determining their 
global behaviour. 

Similarly to the case of the pure YM theory (a = 0) discussed in [^j the global behaviour 
of the solutions can be characterized by the function N{t). Since is positive close to the 
boundary point (r = or r = r/j) it will either stay positive for all r or have a zero for some 
finite tq. In the latter case, which turns out to be the generic one, N has to change sign (a zero 
of even order leads to = 0). From Eq. (|Tl|a) it follows that r has a maximum at tq. Since 

stays negative for r > tq the function r runs back to r = 0, where the solution develops a 
curvature singularity. The corresponding spaces t =const. are compact singular 3-manifolds. 

In the case when N stays positive two possibilities arise. Either r grows without bound or it 
tends to some finite limit. The former case yields the globally (quasi) regular, resp. black hole 
solutions. In the latter case the solution runs into the fixed point A^ = 0,r = l, 1^ = 0. The 
t = const, hypersurfaces are non-compact developing a 'cylindrical throat' for r — > cxd similar to 
the extremal RN solution for r ^ — oo. However in contrast to the latter the metrical function 
A diverges for r — > oo. Thus this solution does not represent the interior of an extremal black 
hole. 

The precise formulation of the described results is the content of Theorem |l]. Its proof, which 
we have deferred to the Appendix, follows closely the one for the analogous Theorem 16 of |^], 
but requires some changes in technical details. 

Theorem 1: For both types of boundary conditions at r = 0, resp. r = the global behaviour 



of solutions of Eqs. (|ljj is completely classified by the following three cases: 

i) N{t) changes sign for some finite value of t, i.e., r(r) attains a maximum. Then N stays 
negative and r turns back to 0. This leads to a compact singular 3-space ('bag of gold'). 

a) N stays positive and tends to 1 for r — > oo and hence r — > oo. In this limit ^ 0, C/ — > 
and K ^ 1. These are the quasi regular, resp. black hole solutions. 

Hi) N stays positive and tends to for r — > oo. The solution approaches the fixed point 
N = U = W = 0,r = K = 1 discussed in the preceding section. There it was shown, 
that for < a < \/3/2 the function W{t) oscillates, whereas for larger values of a the 
behaviour is exponential for r ^ oo. 

The main result based on the above classification is the existence of quasi regular, resp. black 



hole solutions of Eqs. (11). Already from the initial conditions at r = it follows that quasi 
regular solutions can exist only for < a < l/V^. Similarly one obtains a certain bounded 
domain in the (a,r/j) plane, where black hole solutions exist. Although parts of the boundary 
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of this domain are known only through numerical calculations, we are able to establish its 
qualitative features with analytical methods. Similarly to the case a = one obtains discrete 
families of solutions, whose members are distinguished by the number of zeros of the function 
W. 

The precise formulation of our 'Existence Theorem' is 
Theorem 2: 

1. For any < a < 1/V2 and any integer n > 0, there exists at least one global quasi 
regular solution with n zeros of W . For the same values of a there also exists at least 
one 'oscillating solution' with infinitely many zeros ofW and < r < 1 whose asymptotic 
behaviour is given in Eqs. Ijsl 



2. a For any < < 1, a < min(-^, ae(^/i)), and n > 0, there exists at least one global black 
hole solution with n zeros ofW. For a = ae{rh), < rh < 2\/2/3 corresponding black 
holes with a degenerate horizon exist. Under the same assumptions for rh and a there 
exists at least one 'oscillating solution' with infinitely many zeros ofW and < r < 1). 

2-b r/i > 1; For < a < and for any n > there is some rfi^max{o!, n) > I such that black 

hole solutions with n zeros ofW exist for 1 < < T/i,max(cK, n). Furthermore rh^max ~^ cxd 
for a — > 0, i.e., there exist black holes for arbitrarily large if a is taken small enough. 

Remark: For a < \/3/2 one may consider the limit of asymptotically flat (quasi regular or 
black hole) solutions when n tends to cxd. Using the arguments of ^ |l^ one can show that (a 
suitable subsequence of) these solutions converge pointwise to an oscillating solution for r < 1 
and to the extremal RN black hole for r > 1. 

Remark: There is convincing numerical evidence that for a > \/3/2 no black hole solutions 
exist. (Compare Fig. 

The proof of Theorem ^ consists of two main parts. First we demonstrate the existence of 
solutions with a regular horizon, following closely the analogous discussion in j^. We then show 
that the quasi regular solutions as well as the extremal black holes are obtained as limits of 
regular black holes. This approach is necessary since the arguments of cannot be applied 
directly to these cases. It also shows how the regular black hole solution tend to, e.g., quasi 
regular ones as r/^ — > 0. 

The existence proof for regular black holes can be shortly summarized as follows: First we 
partition the initial data into sets corresponding to the three cases introduced in Theorem |l] and 
the number of zeros of W. Using the smooth dependence on initial data we establish that the 
sets Reg„ corresponding to asymptotically flat solutions with n zeros (Case ii) must lie between 
open sets Sing„ or Sing„^^ corresponding to 'singular' solutions (Case i). Analogous results 
hold for solutions corresponding to Case iii. Knowing that the open sets Sing„ or Sing„^^ are 
not empty we conclude that there exist initial data in Reg^. 

We have seen in Sects. |^ and ^ that local solutions with a regular horizon exist as long as 
the condition A^i > is satisfied (compare Eq. (|Tg|)). For any a these initial data form a subset 
of the (r/j, Ty^)-plane shown as region C in Fig. |2|; for given a and there are at most two 
VF/i-intervals. 

Excluding the RN black holes (|l^ with W = and using the invariance of Eqs. (11) under 



the 'reflection' {W,U) {—W,—U) we need only consider initial data with Wh > 0. Based on 
the classification in Theorem |l] and using some results from its proof we can partition the set of 
initial data with A'^i > and Wh > into various subsets as follows. 

For Case i there is some tq such that N{to) = and a corresponding point (ro,Wo) = 
(r(To), VF(to)) in the (r, 14/^)-plane. We introduce Sing^, n = 0, 1, . . . as the set of initial data 
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such the zero of is in region B of Fig. |2| where Wq + o^Tq > 1 and W has n zeros. Sing^^ is 
the set of initial data such that the zero of N is in region A where Wq + a^r^ < 1. 

For Case ii we let Reg„, n = 0, 1, . . . be the set of initial data for regular, asymptotically 
flat black holes with n zeros of W. 

Finally we denote the set of initial data corresponding to Case iii by Osc for a < 
and Exp for a > V3/2. Solut ions with initial data in Osc have infinitely many zeros of W, 
i.e., oscillate. Solutions with initial data in Exp have finitely many zeros and we can further 
partition the set Exp into sets Exp„ according to this number of zeros. As discussed in Sect. 
e'^/'^W is unbounded for generic initial data in Exp when a < 1. There is, however, a subset in 
Exp, corresponding to a lower dimensional stable manifold such that e'^^'^W remains bounded. 

Let us next study the change in the behaviour of solutions induced by a small change in the 
initial data, exploiting the continuous dependence of the solutions on their initial data. 

Case i is generic in the sense that it is stable under sufficiently small changes in the initial 
data. Moreover, for a zero of N in region B the number of zeros of W does not change as long 
as Wh > 0. Therefore Sing„ and Sing^^ are open sets. 

For solutions with initial data in Reg„ the n zeros of W occur for r < 1/a and \W\ mono- 
tonously decreases when r > 1/a. Given any ri ^ the properties of the solution are 
unchanged for r < ri under a sufficiently small change of the initial data. The function \W\ 
may, however, start to increase for r > ri with or without an additional zero of W. Therefore 
a sufficiently small neighbourhood of a point in Reg„ consists of points in Reg„, Sing„, or 
Sing„+i. 

Solutions with initial data in Osc or Exp satisfy > and r < 1 for all r and furthermore 
— > as r ^ oo. Given any ri ^ 1, these inequalities remain valid for r < ri for a sufficiently 

small change of the initial data. They may, however, be violated for some r > ri with either 
= and r < 1 or > and r = 1 while \W\ <C 1. Given initial data in Osc and any positive 

integer n, a sufficiently small change of the initial data yields a solution with at least n zeros of 

W. Therefore a sufficiently small neighbourhood of a point in Osc consists of points in Osc, 

Singoo, or IJ (Reg„ U Sing^). 

m>n 

The analysis for Exp is slightly more complicated. A sufficiently small neighbourhood of a 
generic point in Exp„ consists of points in Exp„, Sing^^, or Sing„. This result is partly based 
on the observation that solutions missing the fixed point have no zeros of W for r > 1, a fact 
directly related to the stability of the RN solution when a > \/3/2 discussed in Sect. |^. The 
neighbourhood of a point in Exp„ for which e'^^'^W is bounded can in addition contain points 
in Reg„, Exp„+i, or Sing„+i. 

In analogy to Lemma 21 and 24 of one can show that for given a and each of the sets 
Reg^ or Osc, resp. Exp consists at most of isolated points. 

Consider a continuous one-parameter family of initial data for some a < 1, e.g., a VF^-interval 
with rh fixed. Prom the above results we can immediately draw the following conclusions: 

1. A family interpolating between Singg^ and SingQ contains at least one point in Osc, resp. 
Exp for a < \/3/2, resp. a > 

2. A family interpolating between Osc and Singg contains at least one point in each Reg^. 

3. A family interpolating between Exp„ or Sing„ and Singg contains at least one point in 
each Reg^ for < m < n. 

We now use the above results for families connecting points near the boundary of the domain 
of initial data with A'^i > and Wh > 0. For a < 1 this boundary consists of three pieces: 

a) The part of the quartic curve A^i = inside the ellipse + ct^^h ~ boundary 
between regions C and A of Fig. ||. 
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b) The part of the quartic curve outside the elhpse, i.e., the boundary between regions C 
and B. 



c) The hue Wh = 0, rh>l. 

For initial data sufficiently close to the curve A^i = we can use the truncated version of 
Eqs. (|T^ ) corresponding to /? = cxd and deduce that there must be a zero of for some r ^ 1 
with r ~ and W ~ Wh- Therefore initial data close to the boundary between regions C and A 
are in Singj^ and those close to the boundary between regions C and B are in Singg. Actually 
all boundary conditions outside the ellipse, i.e., in region C_, are in Singg because Wh 7^ 
implies UW > for all r > 0. This is true in particular for all boundary conditions for a > 1 
except those for the RN solution with W = 0. This concludes the existence proof for regular 
black holes with r/i < 1. 

In order to prove the existence of regular black holes with n zeros of W, i.e., that Reg„ is 
not empty for given a and r/j > 1, it is sufficient to show that solutions with Wh ~ have at 
least n + 1 zeros, i.e., have initial data in IJ (Reg^ U Sing^). Solutions with r/j > 1 remain 

close to the RN black hole (|l4|) on any finite r-interval if Wh is chosen small enough. 

We first want to prove the existence of regular black holes for any r/j > 1 and n. Given 
rh and n we choose a <C For r/j ^ r <C 1/a we may use the approximate form W ~ 

C-y/r sin In r + 9^ valid for 1 k and to argue that the solution has at least 

n + 1 zeros. 

Finally we consider regular black holes for given a < \/3/2 and n. Choosing initial data 
near the curve A'^i = (i.e., max(rft — <C 1) the solution comes close to the fixed point 

r = K = \,N = W = U = ^ and stays there long enough to have at least n + 1 zeros. This can 
be deduced from the approximate form of W given in Eq. (p9[). This concludes the existence 
proof for regular black holes with r/i > 1. 

In the proof of existence for regular black hole solutions we have followed the same strategy 
as in Q which indeed went through with some minor changes in the arguments. Trying to extend 
this stragegy to the quasi regular solutions and extremal black holes we encounter difficulties. 
An essential part of this strategy has been to determine the behaviour of solutions with initial 
data close to the boundary of their allowed region. Whereas for a = the regular solutions 
behave as ~ 1 — near r = 0, this changes to « 1 — + . . . + hr'^ with 7 > 2 for 
the quasi regular solutions with < a < l/\/2 (compare Eq. (17)). The method used in 



to analyze the behaviour of solutions with large h is not directly applicable for a > 0. This is 
particularly evident for a l/\/2 where 7 — > 00. A similar situation is encountered for black 
holes with a degenerate horizon. 

On the other hand it turns out to be possible to obtain the quasi regular solutions and 
extremal black holes as limits of regular black holes. Having established the global existence of 
(asymptotically flat and oscillating) regular black holes, we shall now consider the limiting cases 
corresponding to points on the boundary of the allowed domain in the (a,r/j) plane. 

Guided by the numerical results we shall prove that regular black hole solutions with r/^ <C 1 
consist of two pieces separated by a large r-interval: Starting at the horizon with = the 
solution first approaches the fixed point N = k = Vl — Sa^, W = 1, U = r = 0. Since, however, 
r > the solution necessarily misses this fixed point and eventually runs away along the stable 
manifold for t ^ —00 corresponding to the quasi regular solutions. In the limit r/j — > the two 
pieces get separated by an infinite r-interval and the outer part lies on the stable manifold. The 
same type of argument was used in [^] to show that the Bartnik-McKinnon solutions converge 
to the extremal RN solution for r > 1 as the number of zeros of W tends to 00. 

In order to construct a quasi regular monopole with n zeros of W for some a < l/\/2 we 
take a sequence of regular black holes with n zeros and r/j — > 0. Since the values (Wo, Uo,No, kq) 
at some regular point r = tq are uniformly bounded we can select a convergent subsequence. In 
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order to obtain a convergent sequence of solutions of the autonomous system (|Tl|) we have to 
adjust the variable r. For that purpose we introduce a shifted variable f = t — To{rh) where tq 
is determined by r(ro) = tq. In view of the uniqueness theorem for the initial value problem at 
a regular point the sequence of solutions converges for fixed f to a limiting solution defined for 
all r > 0, representing the outer part mentioned above. 

In order to study the inner part we consider the solutions for fixed r. Introducing p = r/rh 



and T = (W — l)/r we obtain from Eqs. (11,12): 



p = pN , (37a) 
1 

t 



N=l{l-2a'^ - - 2U^ - AT^) + 0{rh) , (37b) 



K = 1 - 20^ - + + 0{rh) , (37c) 
t=U-NT, (37d) 
U = 2T - {k- N)U + 0{rh) . (37e) 

The geometrical structure of the allowed domain of black hole boundary conditions for close 
to zero (see Fig. |l]) implies the existence of a uniform bound for the the initial values = T{0) 
at the horizon. Hence way may select a further subsequence for which the values converge 
to some T, the initial value of the limiting solution. We will show that T must be zero because 
otherwise would vanish for some finite r. In the limit r/^ — > Eqs. (^) imply 

(TC/y = 2T2 + - kTU > (\/8 - k)\TU\ . 

Assuming T 7^ we conclude that TU grows exponentially as soon as k < -v/S, because TU > 
for r > 0. Using the Schwarz inequality it follows from Eq. (|37|b) that must vanish for some 
large enough r, contradicting the properties of the limiting solution. Therefore T = and hence 
T = U = 0. Solving the remaining equations we conclude that the limiting solution tends to the 
fixed point N = k = Vl — 20^, = 1, C/ = r = for r ^ 00. 

The remaining argument closely follows the one used in |jl^. The solutions for sufficiently 
small rfi come arbitrarily close to the fixed point, miss it, and run away along the stable (for 
r —00) manifold. In a small but fixed neighbourhood of the fixed point the distance from 
the stable manifold decreases exponentially with increasing r. This implies that in this small 
neighbourhood the outer part of the limiting solution lies on the stable manifold, i.e., is a quasi 
regular solution. 

The arguments for the existence of extremal black holes are fairly similar to those for quasi 
regular solutions described above. The values of r, W, and k at a degenerate horizon for 
l/\/2 < a < 1 are 



V2a2 - 1 ^ 



We = — — — , He = V2a^ - 1 , (38) 
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and {Wh — We)l \jTh — re is uniformly bounded. This suggests to introduce new variables p = 
{r -re)/{rh-re), n = N/{r-re), w = {W -We)/ y^re{r - re), and u = rV j ^Jveir - re). Given 
a sequence of (asymptotically flat or oscillating) non-degenerate black holes with — > re we 
can select a subsequence such that {Wq, Uo,No,ko) at some regular point rg > re as well as Wh 



converge. Rewriting Eqs. (11) in terms of the rescaled variables p, n, w, and u and neglecting 



contributions that vanish as r/i — > rg we obtain 



p = pn , (39a) 
h= {k — n)n — 2u^ , (39b) 
k = Ke — , (39c) 
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w = —-nw + u, (39d) 

ii = 2W^w -{k + ^n)u , (39e) 

with the constraint 

Kn = Kl + u^ - 2W^w'^ . (40) 

Solving Eqs. (pg|) with initial data p{0) = 1 and w{0) = Wh, we first observe that k ^ k^, and n 
remains bounded as r ^ oo. Next we want to show that w and u remain bounded. Note that 
wu > for r > except when Wh = and hence w = = u. For large r when hk, the linear 
system for yTpw and y/pu has the two eigenvalues 

71,2 «^(-«e± 7/^2 + 8^2)^ (41) 

and therefore y^ti^ and -^pn both grow as e^'^'^ . This is, however, not possible if 271 > Kg since 
(pn)' = Kpn — 2pv? would then imply that n vanishes for some finite r. Therefore w and u 
are either bounded (for 271 < Kg) or are identically zero (for 271 > Kg). We conclude that the 
inner part of the solution for — > approaches the fixed point r = rg, W = We, k = Atg, and 
N = U = Oasr— >oo. The remaining step of the argument proceeds as above for the quasi 
regular solutions. 



6 Numerical Results 

The problem to find global solutions of Eqs. (P) by numerical integration is complicated by the 
singular nature of the boundary points. As was discussed in Section Q the regular solutions lie 
on submanifolds of the phase space, the 'stable manifolds' of these singular points. Hence the 
construction of global solutions may be reduced to finding intersections of the stable manifolds 
emanating from the two boundary points. This, however, requires to extend the stable manifolds 
from the neighbourhood of the boundary points to some common value of r. A general strategy 
to achieve this goal is the following. Close to the singular points one replaces the differential 
equations by a system of integral equations implementing the correct boundary conditions at the 
singular point ||T^. The stable manifold is then parametrized by a subset of the boundary values 
at the other endpoint of integration. The integral equations, resp. suitable discretizations, may 
be solved by iteration (a procedure guaranteed to converge close enough to the singular points). 
The output of this first step is a set of initial data for the differential equations at regular points, 
which may then be integrated using, e.g., the Runge-Kutta method to some common value of r. 
Varying the free boundary values of the integral equations at either end one may try to intersect 
the two stable manifolds at this common value of r. 

The procedure just described is admittedly rather complicated, yet certain simplifications are 
possible in most cases. Since the numerical integration procedures come with certain numerical 
errors anyway, one may as well approximate the stable manifolds close to the boundary points 
by their linearizations, provided this does not introduce unduely big errors. This allows to 
circumvent the use of the integral equation. Another case in which the integral equation can be 
avoided is, that the stable manifold can be directly parametrized at the singular point. In this 
case one may use the Runge-Kutta method immediately from the boundary point. Thus in the 
simplest case one starts, e.g., directly at r = 0, resp. r = integrating the Eqs. (^) with the 
help of the Runge-Kutta method as close as possible to the other boundary point, e.g., r = 00. 
Linearization close to this singular point results in unstable modes depending on the choice of 
the free boundary values at the starting point. Studying the effect of finite or infinitesimal 
variations of these values one can try to suppress iteratively the unwanted unstable modes, thus 
improving the solution. Cases where this simplified approach fails are extremal black holes, the 
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Figure 4: The fundamental and first excited monopole solution for {3 = oo and a = l/\/2 
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Figure 5: The parameter h vs. a of the fundamental monopole, the first two excited ones, and 
the oscillating, resp. exponential solution for /3 = and /9 = 4 



'quasi regular' solution for a close to l/\/2 and solutions for large but finite values of /3. The 
reason for this failure in the last case is the rapid growth of the unstable Higgs-mode (~ e^^) 
before one reaches the asymptotic region for W . 

Fig. ^ shows the fundamental and first excited quasi regular solution (/3 = oo) for the 
extreme case a = 1/ \/2. These solutions could not be obtained without using integral equations 
to determine the stable manifold at the origin. 

Let us now proceed to a qualitative description of the domain in the space of the parameters 
Q, /?, r/j and Wh for which global black hole solutions exist. Starting from one particular 
solution one may vary these parameters and try to find a whole family of such solutions in a 
connected part of the parameter space. In view of the smooth dependence of the solutions on the 
parameters the boundary of the corresponding region in parameter space must be determined 
by some kind of singular behaviour. A typical phenomenon of this type is that the solutions on 
their way from r = r/j to r = oo get caught by some fixed point different from r = oo. This is 
indicated by the fact that the solutions spend more and more 'time' (in r) in the vicinity of this 
wrong fixed point. A second phenomenon is a change in the nature of the boundary condition 
at r/j. A trivial boundary of this type is given by the condition > 0. Another case is the 
degeneracy of the horizon. 

Beginning with the case /? = already described in some detail in I, it is instructive to 
observe what happens with the fundamental monopole, if one starts from the flat space Prasad- 
Sommerfield solution and increases a (in order to obtain the flat limit a rescaling r ^ r/a is 
required and hence 6(a) ~ for small a). As seen from Fig. ^ the parameter b increases 
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Figure 6: Domains of existence £ for non-abelian black holes: a) for /3 = 0, 1, 2, 3, and 4; b) 
for /? = 6 and cxd 



monotonously while (y. runs through, a maximal valuG o^max* 

Finally the regular monopole ceases 
to exist for some critical value of a = Qc- At this point the solution runs into the 'wrong' fixed 
point with = at r = 1, i.e., bifurcates with a solution belonging to Exp. For a between Oc 
and ttmax we have two different monopole solutions with a bifurcation point at Omax- As usual 
this bifurcation is accompanied by a change in stability. Whereas the solutions on the lower 
branch of the curve h{a) are stable against linear perturbations the ones on the upper one have 
exactly one unstable mode . It is suggestive to relate this unstable mode to the gravitational 
instability observed for the Bartnik-McKinnon (BM) solution |17|. In fact, as seen from Fig. ^ 
there is a second branch of the h{a) curve for the first excited solution starting with the BM 
value h ~ 0.4537 for a = and ending at the critical value Uc = \/3/2. All the solutions on this 
branch have also one unstable mode. Similar curves exist for the higher excited solutions, all 
of them ending at Oc = \/3/2 with the same value of 6(ac). At this common limiting point all 
the excited monopole solutions bifurcate with the solution in Osc, i.e., run into the fixed point 
with AT = at r = 1. 

The same type of phenomena persists for black holes of small radius r^. For given non- 
abelian black holes exist only for a finite a-interval. The domain of existence £ in the (a, r/j) 
plane is shown in Fig. ^ (where avh is plotted versus a). The difference between amax and Uc 
for the fundamental solution decreases with increasing and disappears for some limiting value 
Th ~ 0.793. Also the value of Oc decreases with increasing tending to Oc = \/3/2 for « 0.991 
still below the maximal value of = 1. This means that the curve delimiting the domain 8 
has a short straight section between the points (\/3/2, 0.991) and (\/3/2, 1). For \/3/2 < a < 1 
the points on the curve {oic{rh)^fh) are determined by the condition that the solution in Exp 
ceases to have a zero as r/j is increased from zero. The reason is that at the critical value of a 
the monopole solution bifurcates with the solution in Exp with e^^'^W bounded (compare the 
discussion in the preceding section). For a = \/3/2 this just happens at Vh ~ 0.991. 

As was already explained in I the curve delimiting the existence domain 8 in the direction 
for a < \/3/2 is determined by a different mechanism. As is varied (for some fixed value of 
a) the value for decreases from 1 to 0. It turns out that the curve {r^, Wh) can be continued 
smoothly to negative values by the reflection Wh — > — W/j, i.e., runs through an extremum 
as Wh is varied from positive to negative values. For small values of a the extremum of rh at 
Wh = turns out to be a minimum, i.e., the maximal value of occurs for some \Wh\ > 
(compare Fig. 0). With increasing a the extremal points of r/, approach each other and merge 
for some particular value of a ~ 0.732. Qualitatively the same phenomenon repeats itself for 
all the excited solutions, the extremal point of Vh with Wh = moving inwards towards rh = 1 
with increasing order of the excitation. 
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For Wh = the non-abelian black holes bifurcate with the (non-extremal) abelian RN black 
hole. As a tends to \/3/2 the corresponding RN black hole becomes extremal. All the extrema 
of r/i correspond to bifurcation points, accompanied by a change of stability of the solution. 

When P is increased from zero the basic phenomena persist and the domain £ changes in a 
smooth way with some clear tendencies. The value of for which the difference between ttmax 
and Oc vanishes moves down quickly from its value ~ 0.991 for /3 = to zero for /3 ~ 0.757. 
Hence for (3 > 0.757 the phenomenon amax 7^ Oc disappears. The same is true on the other 
limiting curve of the domain £. The value of a where the extrema of rh merge decreases with 
increasing /3 and becomes zero for /? « 1.18 for the fundamental monopole (and similarly for the 
excited solutions). The values of r/j with Wh = 0, however, turn out to be /3 independent, i.e., 
for P > 1.18 the boundary curve of the domain £ for a < \/3/2 is independent of p. Fig. |^ 
shows the domains for some values of /5 < 4. 

A new phenomenon appears for /? > 4, because the boundary of the allowed domain of 
the (a,r/i) plane for black hole boundary conditions now intersects the region a < \/3/2. The 
curve delimiting the domain £ towards growing a develops further branches as compared to 
/? < 4. A picture for /3 = 6 is shown in Fig. ^3. Starting from the bifurcation point with the 
extremal RN solution at a = \/3/2, r/j = 1 the boundary curve is the straight line a = \/3/2 to 
the intersection with the curve for an extremal horizon. Then it follows the latter curve to its 
extremal point P with the minimal value of a. The continuation of the curve from there with 
decreasing and increasing a is determined by a new phenomenon. The solution runs into the 
fixed point with = corresponding to the conditions (pO|). For values of /3 somewhat bigger 
than 4 this branch of the boundary curve is valid up to a point Q with a = V3/2. The value 
of r^i at Q and the segment PQ in Fig. ^3 are drawn schematically due to lack of numerical 
data. From Q there follows again a straight section with a = \/3/2 till the curve adrh) already 
discussed for /3 < 4 branches off. The latter is then relevant down to r/i = 0. With increasing /3 
the whole boundary curve moves inwards approaching the limiting curve for /3 = oo also shown 
in Fig. ^b. 

7 Stability of RN 

It is well known |l8[ that the Reissner-Nordstr0m black hole is a stable solution of the Einstein- 
Maxwell theory. Allowing for the more general perturbations provided by a non-abelian gauge 
field the question of stability has to be reconsidered. It was shown in |l^ that the RN black 
holes, considered as solutions of an EYM theory, are unstable. The situation is slightly more 
complicated when a Higgs field is included. The Higgs field itself does not contribute unstable 
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modes, moreover it has a stabilizing influence through the mass of the Yang-Mills field due to 
the Higgs eff'ect. It is, in fact, easy to see that for values of the YM mass a > 1/r/i no instability 



of the type considered in |19] occurs. The precise value a{rh) beyond which the RN solution 
becomes stable can, however, in general be determined only numerically ||3|, except for the 
extremal RN black hole with = 1. We will show that in the latter case the solution changes 
stability exactly for a = \/3/2. More precisely we will prove that for a < \/3/2 there is an 
infinite number of unstable modes, all of them disappearing for a = \/3/2. This is in accordance 
with the numerically observed bifurcation of the extremal RN black hole with the non-abelian 



ones (for all n) and the known change of stability at bifurcation points [20|. 

In order to investigate the stability we consider the linearized field equations in the RN 
background. It is easy to see 110], that the metric perturbations vanish and it suffices to study 



the linearized equation for W. With the ansatz W = e w this eq. reads as 



+ (aV2 - l)w + (1 - -)— = io^- -^w , (42) 



dr^ r dr (r — 1) 

with r = 1 + e^. Introducing a new coordinate p with 



„2 



dp = -dr , (43) 

(r-1) 



Eq. ( p2[ ) becomes a standard Schrodinger equation 

d'^W 9 1^9,9 1 



,^+Vw = u;^w with y = (1 - -)^(a^ - — ) . (44) 
dp^ r 

Note that the coordinate p tends to — oo for r — > 1 and to +00 for r — > 00. Near the horizon 
r = 1 the potential V behaves like {p? — l)//0^, hence assuming < 3/4 we get V < —l/Ap^ for 
p — > —00. According to a standard text-book theorem p.l463ff, this implies the existence of 
infinitely many bound states for the Schrodinger eq. ( ^^ accumulating at uj"^ = 0. This proves 
the instability of the RN solution for a < \/3/2. 

For the opposite case > 3/4 we shall now show that there are no bound states using the 
so-called Jacobi criterion. The corresponding condition is, that the zero-energy wave function w 



with the appropriate boundary condition at r = 1 (guaranteeing regularity) has no zero |21, 22]. 
Actually it suffices to study the limiting case = 3/4 in view of the monotonicity of the 
spectrum of the Schrodinger operator with V. 

It turns out to be convenient to use the coordinate R = r — 1 instead of p in this case. 
Putting V = ^/Rw/r the zero-energy Schrodinger equation reads 

R{Rv'y = 4(/+i)2 (-^' + + 5^ - \> (45) 
We transform it into a Riccati equation putting y = Rv' jv 



The shift y ^ z = y^ R/2 leads to 



(47) 



2(i? + l)2 R 



Regularity of the solution at i? = requires y(0) = and hence also z(0) = 0. From Eq. (|47| ) it 
follows that z > for all ii > and thus y > —R/2. For large R the r.h.s. of Eq. ( ^6[ ) is positive 
as long as —R/2 < y < and hence y is bounded from below. This shows that v and thus w 
have no zero, implying stability. 
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Appendix: Proof of Theorem |I] 



Before entering into the proof of Theorem |T| we make some general observations concerning the 
behaviour of the solutions. From Eq. (pT|e) it follows that W{t) can have no maxima (minima) 
for W > {W < 0) outside the ellipse W'^ + a^r^ = 1 in the (r, W) plane. Hence a solution 
with UW > outside this ellipse cannot turn back to W = for increasing r. 
Using the constraint Eq. (^) we may rewrite the Eq. (pA[b) for in the form 



iV = -(l-iV2-2C/2-^ _L^_2a2p72N /4g^ 

2 



from which we see that if < 1 holds for some r it remains true for all larger r. Combining 
Eq. (^8|) with the one for k we get 

(k- N)-= -k{k- N) + 2U^ + ^ ' (49) 

It follows that if K > for some r it holds for all larger r. Introducing rj = n + N — 2 we find 

f] = -r] - ^7]'^ - ^{k - Nf - 2aV , (50) 

implying k + A^ < 2 for large enough r. 

Another useful inequality follows from the equation 

(1 - kNJ = -N{1 - kN) + 2{k - N)U^ + 2a^NW'^ . (51) 

Since kN < 1 at the origin, resp. horizon this inequality holds for all r. Together with the 
constraint Eq. ( [l^ ) this implies the inequality 

2C/^ < 1 + ^ ^3-^ + 20^14/^ . (52) 

Hence U stays bounded for r > as long as W is bounded. 

Let us now come to the proof of Theorem [l|. It runs essentially along the same lines as 
the one given for a = in therefore we only indicate the necessary modifications of the 
arguments given there. 

For Case i we observe that the constraint Eq. (IT^ ) implies the relation 
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+ 2a^W^ = 1 + 2C/o^ > 1 , (53) 



where Wq, Uq and vq denote the values of W, U and r at the position of the zero of A^. The 
domain of pairs (ro, Wq) allowed by this inequality is depicted as regions A and B in Fig. |2| (for 
a = 0.75). From the preceding discussion it follows that this domain is just the complement of 
the one allowed for black hole boundary values. 

Whenever W'^ — 1 + a^r"^ > and UW > for some r one necessarily ends up with Case i. 
The proof of this statement is completely analogous to the one for a = given in Ref. In 
view of the assumed boundary conditions, when a > 1 one obtains Case i (except for the RN 
solution with W = 0). Therefore, in order to avoid that N{t) develops a zero for finite r, the 
black hole boundary conditions have to be taken within the intersection of the domain C and 
the interior of the ellipse {r^, Wh) : + a'^rf^ = 1. This region will be denoted by C+. 

The argument showing that r(r) turns back to r = within a finite r- interval, if N{t) has 
a zero, can be taken more or less literally from Q. 
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Now let us turn to Cases ii and iii. Without restriction we may assume a < 1. In these cases 
N(t) stays positive and r(r) grows monotonously, hence r(r) is either unbounded or tends to 
some finite limit for r — > oo. 

In Case ii r(T) is unbounded and we may assume that UW < for large r (otherwise we 
end up with Case i). From Eqs. ( [ll] ) we get, however, 

— 1 

{UWy = rU^ + + aV) - {k - N)UW , (54) 

r 

implying (UW)' > (remember k — > 0). Taking into account VF^ > and {W^J = rUW, 
we conclude UW — > for r oo. This in turn implies U,W ^ and k,N ^ 1. 

There remains only to discuss Case iii when r(r) stays bounded. As argued in Q this 
implies iV — > for r — > oo. From the contraint (|l^) it then follows that the solution has to stay 
in the closure of region A in the limit as r — > oo. Since 1 — 2o?W'^ has a positive lower bound 
in region A (except when a = l/\/2 and r = 0) one obtains a positive lower bound for k when 
T is sufficiently large (see Eq. (pT|c)). 

To prove that the solution runs into the fixpoint U = W = as claimed, we shall employ 
the following Ljapunov function (for the {U, W) subsystem) 

h = ^W^ + ^W^{l-a^r^ -^W^) + 6WW . (55) 

We may consider /i(r) = + (1 — a^r^ — VF^/2)y^/2 + 6xy as a quadratic form in {x,y) = 
{W, W), being positive in or close to region A for < 6 < a/(1 — a-^) /2 and vanishing only for 
W = W = 0. The derivative of h is given by 

h = -((1 - aV^ - W'^)6 + a'^r^N)W^ - {k - 2N - 6)W'^ - (k - 2N)6WW , (56) 

which is a negative definite quadratic form for large r provided 5 is chosen small enough. Since 
two definite quadratic forms are relatively bounded there exists some constant c > such that 
h < —ch for large enough r. Hence /i(t) — > and thus also C/, — > for r oo. From 
Eqs. (11c,|12|) one deduces that then k and r tend to 1, as claimed. This concludes the proof of 



Theorem |lj. 
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